Chaotic dynamics in superconducting nanocircuits 
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The quantum kicked rotator can be realized in a periodically driven superconducting nanocircuit. 
A study of the fidelity allows the experimental investigation of exponential instability of quantum 
motion inside the Ehrenfest time scale, chaotic diffusion and quantum dynamical localization. The 
role of noise and the experimental setup to measure the fidelity is discussed as well. 
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The kicked rotator is a paradigm model in classical 
and quantum chaos. The chaotic regime of the classi- 
cal rotator is characterized by the exponential separa- 
tion of nearby trajectories, with rate given by the maxi- 
mum Lyapunov exponent A, and chaotic diffusion in the 
momentum (action) variable P|. The quantum kicked 
rotator (QKR) is exponentially unstable only up to the 
Ehrenfest time Ie needed for a minimal quantum wave 
packet to spread in the angle coordinate of the action- 
angle phase space • The classical-like diffusive behavior 
is destroyed by quantum interference effects, leading to 
a dynamically localized state after the localization time 
^* 0,3- Since typically t* ^ Ie, the diffusive behavior is 
possible also in the absence of exponential instability. A 
quantitative description of classical to quantum crossover 
have been recently presented 0|. 

Despite the long-standing interest in the QKR, only 
few proposals have been put forward and the only ex- 
perimental implementation so far has been realized with 
cold atoms exposed to time-dependent standing waves of 
light 0. In this Letter we suggest, for the first time, 
a way to realize the QKR by means of a superconduct- 
ing nanocircuit 0- We analyze the chaotic dynamics 
of a periodically driven Superconducting Single Electron 
Transistor (SSET) and show that under appropriate con- 
ditions it reduces to a "generalized" QKR as the external 
phases of the superconducting electrodes can be used to 
tune the quantum dynamics of the superconducting de- 
vice. 

A way to quantify the stability of quantum motion 
is to study the fidelity it measures the overlap of 
two states obtained through two slightly different evolu- 
tions 0, [11 [II El H [li. Following the idea put for- 
ward by Gardiner et al. [l5j we also discuss how it is pos- 
sible to measure the fidelity for our proposed Josephson- 
QKR. The flexibility in the design of superconducting 
nanocircuits allows us to consider several different situa- 
tions. In the semiclassical regime and for strong enough 
perturbations the fidelity decay, exponential with rate 
given by the Lyapunov exponent or power-law, follows 
the classical one up to the localization time t*. For t > t*, 
the fidelity oscillates around a value given by the inverse 



of the localization length. 

The system we consider, illustrated in Fig. [H is very 
closely related to the Cooper pair shuttle [la, UM but it 
operates in the regime where the Josephson coupling is 
much larger than the charging energy. As it will be dis- 
cussed later, the capacitive coupling to a Cooper pair box 
is needed for the measurement of the fidelity. The Cooper 
pair shuttle is a superconducting device composed by 
a small superconducting island coupled to two macro- 
scopic leads The couplings to left (L) and right (R) 
electrodes are time dependent with period 2T and the 
island is never connected to both leads simultaneously. 
The two leads are macroscopic and have definite phases 
4'L,R^ while the superconducting island is described by 
the number n of excess Cooper pairs present on it. In 
the Cooper pair box only states consisting of zero (|0)) 
and one (|1)) Cooper pairs are allowed. The Hamiltonian 
describing the system is 
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where if is the conjugate phase of n {[n,ip] 
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Ec = (2e)2Ca/(CsCa - t) is the charging energy, Cs 
and Co- are the total capacitance of the Cooper pair shut- 
tle, and of the Cooper pair box, the dimensionless gate 
charges are defined as rig — VgCg/{2e), nc = VcCc/i^e). 
The condition C < Co- <C Ce guarantees that in the 
Cooper pair box, only states |0) and |1) are relevant. We 
set to zero the Josephson energy of the Cooper pair box 
and < fj, = Cint/Ca < 1. The time dependence of the 
Josephson energies Ef'-^\t) are plotted in Fig.[T] When 
the island is coupled to one of the leads ("Josephson 
kick" ) the corresponding Josephson coupling has value 

Ej, otherwise E^f\t) = Ej^\t) ~ 0. We employ a sud- 
den approximation (switching time At <C ^/Ej) so that 

E^j^'^\t) can be approximated to step functions. 

We first set /i = and study the chaotic dynamics of 
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Figure 2: Upper figure: {(An)^) = ((n — {ii))^) as a function 
of time for K = 10, k = 15, /i = 0, and phase difference (from 
bottom to top) (/> = 0,0.05,0.1,0.4,0.8,0.25 in units of 27r. 
Lower figure: Saturation value of {(An)^) oc i'^ as a function 
of <j} for tfie same parameter values as in tfie upper figure. 



Figure 1: Upper panel: Schematic drawing of a Cooper pair 
shuttle (dashed red box) capacitively coupled to a Cooper pair 
box (dot-dashed green box). Lower panel: Time dependence 
of the left and right Josephson energies within a single period 
2T = 2{tj + tc). 



the Cooper pair shuttle (dashed red box in Fig. 01. For 
the sake of simplicity we assume Vg = Vq = 0. The 
Hamiltonian reduces to that of a QKR 

Ho = -EjY,[cos{'p-MSit-2nT) 



+ cos{ip^ <j)B)6{t- {2n+l)T)] 



(4) 



if the effect of charging energy during the "Josephson 
kick" can be neglected. This condition, under the as- 
sumption that Ej 3> Ec, is satisfied if the charging 
term cannot induce a significant change of ip during 
the kick |l8l |. This sets a limit on the maximum num- 
ber of allowed charge states involved in the dynamics 
{nEc'tj/h < 1). The dynamics of the Cooper pair shut- 
tle mimics that of a QKR with the additional free param- 
eter (f) = 4'R~ 4'L [13. As the parameters k = Ejtj/h 
and K = {Ectc /K){Ejtj /K) are varied, the dynamics 
of the QKR exhibits several interesting phenomena, in- 
cluding quantum ergodicity, quantum resonances and dy- 
namical localization 0. The classical limit corresponds 
to fc — > 00, with K = const. The classical dynamics 
corresponding to Eq. I^J depends only on the parameters 
K and (f). For A' > 1, the dynamics of the charge on 
the central island is diffusive: ((77.2* — no)^) *— ^ D{2t), 
where D is the diffusion coefficient and t is time measured 



in units of T as it will be in the following. Following 
Rei. Enjl ■ we obtain 



D = - 
2 



1 - 2cos(20)J2(i^) 



(5) 



In the semiclassical regime k ^ 1, the QKR follows the 
classical diffusive behavior up to the localization time 
t*. For t > t* quantum interference effects, as shown 
in Fig. 13 (upper panel) suppress this chaotic diffusion: 
The wave function is exponentially localized in the charge 
basis, over a localization length £ (f ~ i* ~ D) The 
fluctuations of the charge in the central island saturate. 
The localization length can be further tuned by changing 
the phase difference as shown in the lower panel of Fig. |2] 
We now turn to the discussion of the fidelity de- 
fined as Fit) = \f{t)\^ = |(7^(0)| /(t) |V(0)) P, where 

/ ~ exp (^j-Hitj cxp (^—jrHQt^ is the so-called echo op- 
erator {Hq and Hi are defined in Eqs.(|21 El)- Under 
the action of this perturbation the fidelity has a sim- 
ple interpretation: The state evolves with the unper- 
turbed Hamiltonian Hq for a time t, it is shifted by 
an amount ^ along the coordinate n and evolves back- 
wards in time with the same Hamiltonian Hq for the 
same time t. F{t) measures the overlap of the final 
and initial state. Note also that the specific form of 
the perturbation implies that f{t) is a 27r-periodic func- 
tion of fiK/k. Although the fidelity depends on the spe- 
cific perturbation Hi — J?n, its time dep endence shows 
rather general features HE O El El Q- At very 
short times t < t^, the fidelity is quadratic with time 
1 — F{t) oc (/it)^ [Sj. At later times the fidelity decays 




Figure 3; Fidelity dependence on time for K — 10, k = 15, 
<j> = 0, — 0.5 and different noise strengths: 'ytc ~ 0, 10"'' 
(upper and lower black curves), 10~^ (green), 10~^ (red). 
Circles and triangles show the decay of fidelity in the clas- 
sical limit (with no noise) and in the limit of strong noise 
respectively. Dashed and full blue lines show the saturation 
of fidelity in the localized regime {F{t) = l/£) and the power 
law decay F{t) oc l/VSt'. 



until it reaches a saturation value l/£ for t > t*, being 
£ the localization length, i.e. the total number of quan- 
tum levels involved in the dynamics. This is shown in 
Fig. 131 where the saturation value provides an indirect 
measurement of the localization length. For tp < t < t* 
general arguments 0, 0| allow to define a perturbative 
border such that for fi < fip 1/V£ the fidelity decay 
is Gaussian. For < u < fic ^ ^ the decay is expo- 
nential with rate ^ /i^ [lO|. For fi > fic, the quantum 
fidelity decays as in the classical case up to the localiza- 
tion time scale t* In particular for times smaller 
than the Ehrenfest time, the decayis exponential, with a 
perturbation-independent rate A |9j, where A is the Lya- 
punov exponent characterizing the exponential instabil- 
ity of classical chaotic dynamics. This behavior can be 
observed in Fig. 0] up to time ~ (1/A) ln(fc//iA') and it 
is followed by a square root decay: F{t) oc l/\/Dt (see 
Fig. El. As can be observed in Fig. |21and Fig. ^ in our 
system the fidelity decay follows the classical one up to 
the localization time scale for experimentally accessible 
values of /i 0.8. 

All the features of the fidelity discussed so far can 
be measured by adapting the protocol presented in 0| 
as sketched in Fig^ By preparing the system in the 
initial state |V'(0)) ® (|0) -I- |1)) /V2 and by applying a 
7r/2-pulse to the Cooper pair box at time t, the fi- 
delity can be extracted by measuring the probability 
Pi = [1 — 3f?e [/(t)]] /2 of the Cooper pair box being 
in the state |1). By repeating the procedure for the 
initial state |'0(O)) (g) (|0) + i |1)) /V2, one can measure 
Pi = [1 — 3m [/(t)]] /2 and therefore the fidelity ampli- 
tude. 




Figure 4: Fidelity decay for K ^ 10, k ^ 2.6 x 10^, 4>l = 
(j>R = 2tt/3. Different curves correspond to different values 
of of /i and "ftc, in the Fermi golden rule regime (diamonds, 
from top to bottom fi = 5x 10"^, 10"\ 3 x 10"\ 5 x 10"^) and 
in the Lyapunov decay Regime (circles, fi = 0.8, 1, 1.5). Black 
triangles are obtained in the presence of noise (ytc ~ 1), for 
fj. = 5 X 10"^3 X 10~'\5 X 10~\ 1. The red line shows the 
exponential decay f{t) — exp(— At) where A is the Lyapunov 
exponent of the kicked rotator for K — 10: \ — 1.62 ~ 
ln(A'/2). 



An important issue to consider is the effect of the ex- 
ternal environment on our system. We focus on effect 
of noise due to gate voltage fluctuations. It amounts in 
adding a term to the Hamiltonian of the form 



(6) 



^(t) being Gaussian distribute with (C(O)stoc ^ ^^"^ 
(C(*)C(*')>stoc = fi^l^it - t')/T. Due to the condition 
Ej ^ Ec, this type of noise is relevant only between 
kicks [U. At the classical level the presence of ^(t) does 
not significantly affect the classical diffusion coefficient 
(see Fig. O. If noise is weak (jtc ^ 1), it gives small 
corrections to the results discussed so far. The effect 
of noise on the fidelity are shown in Fig. O and Fig0] 
resulting in the destruction of the dynamical localization. 
In the limit of jtc ^ 1, the system reaches the classical 
behavior characterized by an exponential decay of fidelity 
with the Lyapunov exponent at short time scales (Fig.QJ 
and by a behavior F{t) cx l/^J Dt at arbitrary long time 
(Fig.©. 

Further insight in the effect of noise can be obtained by 
means of master-equation approach which in some lim- 
iting case allows for an analytical treatment. The noise 
introduced in Eq. @ can be traced out and the evolution 
of the reduced density matrix is given by: 



miT 



7 /.2- 



2nph + ph 



(7) 



The diagonal terms of p, pn = {n\ p\n), are not modified 
by the presence of noise during the free evolution time. 
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while the ofF-diagonal ones are exponentially suppressed 
on a time scale ^7"^. Therefore, as 7 ^ cx), only diag- 
onal elements p„ of p survive, and they are determined, 
at integers multiples of periods, by the map 

Pm(i+l)=^J^_„(fc)pn(<)- (8) 
n 

The map (0 gives a diffusive behavior: {{n^ {i))) ^^^^ = 
= (fc^/2)t. This justifies the classical decay 
of the fidelity, F{t) oc l/^/Dt, discussed above. Note 
also that in the strong damping limit the diffusion coef- 
ficient is independent of (j). Noise destroys time correla- 
tions (sin(0t) sin(0t-i-2)), from which the dependence on 
arises. 

The set-up proposed in can only access the av- 
eraged (over the noise) fidelity amplitude. Except for 
the case of weak noise, the averaged fidelity amplitude 
displays rather different behaviors with respect to the 
averaged fidelity. The evolution of the fidelity ampli- 
tude can be determined by observing that f{t) ~ Tr[f] 
and that / fulfills the same differential equation as the 
density matrix, Eq. (O, once the replacement [H^p] — > 

(^fHo — Hif^ has been performed. The same argument 
used for the density matrix leads to the conclusion that 
in the 7 00 limit / is diagonal and its evolution is 
described by the map 

f„,{t + 1) = e-*f Jn.-nik)Mt) , (9) 

n 

where /„ = (n|/|ri). Using this map, the asymp- 
totic decay of |(/(i))stocl '^^^ ^'^ computed analyti- 
cally and we obtain \{f{t))stoc\ ~ exp(— ci), with c = 
-^/q ln|Jo[2A:sin(6»/2)|d6'. In measuring the fidelity 
amplitude also noise effects due to fluctuations of the 
Cooper pair box's gate voltage have to be taken into 
account. These last fluctuations are uncorrelated to 
the previous one, and are treated in the same way by 
adding the term |f) (f | to the Hamiltonian. One gets 

(^.>stoc + ^iP9')stoc = (1 + - eM-mfit))stocP, 

where F is defined through (5(<)S(t'))stoc = fi^T6{t - 
t')/T^. 

Finally we would like to comment on the experimental 
feasibility of our proposal. Due to physical constraints, in 
the proposed setup we cannot explore the whole parame- 
ter space. By choosing tj ^ fO~"'^°sec and Ec ^ f 0~^cV, 
we can access parameter values corresponding to interest- 
ing physical regimes. For instance, as shown in Fig.|21we 
can observe dynamical localization for K = fO, fc = f 5 
(corresponding to ~ 5 x 10~^sec, Ej ^ fO^'^eV). The 
semiclassical regime and the Lyapunov decay can be ob- 
served for K = 10, k = 2.6 x 10"^ (sec Fig.O correspond- 
ing to ~ X 10~^sec and Ej ^ 10~^eV. For this choice 
of parameters the maximum number of levels for which 



the QKR correctly describes the physics of the system is 
fi/{Ect.j) ~ 6 X 10^. 
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